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Abstract. In this paper we describe a new data structure that supports orthogonal range reporting 
queries on a set of points that move along linear trajectories on a U x U grid. The assumption that 
points lie on a U x U grid enables us to significantly decrease the query time in comparison to the 
standard kinetic model. Our data structure answers queries in 0( \/log U / log log U + k) time, where 
k denotes the number of points in the answer. The above improves over the [2(lagn) lower bound 
that is valid in the infinite-precision kinetic model. The methods used in this paper could be also of 
XJ} . independent interest. 
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O \ 1 Introduction 

Data structures for querying moving objects were extensively investigated in computational geome- 
try and database communities. The orthogonal range reporting problem, i.e. the problem of storing 
a set of points S in a data structure so that all points in a query rectangle Q can be reported, was 
also extensively studied for the case of moving points. In this paper we describe a data structure 
that supports range reporting queries for a set of moving points on a U x U grid, i.e., when all 
, point coordinates are positive integers bounded by a parameter U. 

Previous and Related Results. The kinetic data structure framework proposed by Basch et 
al. [7] is the standard model for studying moving objects in computational geometry. The main idea 
of their approach is to update the data structure for a set S of continuously moving objects only 
at certain moments of time: updates are performed only when certain events changing the relevant 
combinatorial structure of the set S occur. For instance, the data structure may be updated when 
the order of projections of points on the x-axis changes or the closest pair of points in S changes; 
see e.g., |7|15j for a more detailed description. 

The kinetic variant of the range tree data structure was presented by Basch, Guibas, and 
Zhang (Sj; their data structure uses 0(n log d_1 n) space, answers d-dimensional queries in 0(\og d n+ 
k) time, and can be updated after each event in 0(log d n) time; henceforth k denotes the number 
of points in the answer. The two-dimensional data structure of Agarwal, Arge, and Erickson [2] 
supports range reporting queries in 0(log n+k) time and uses 0(n log nj log log n) space; the cost of 
updating their data structure after each event is 0(log 2 n). As follows from standard information- 
theoretic arguments, the O(logn) query time is optimal in the infinite-precision kinetic model. 
Linear space kinetic data structures were considered by Agarwal, Gao, and Guibas [3] and Abam, 
de Berg, and Speckmann [lj. However these data structures have significantly higher query times: 
the fastest linear space construction [1] answers cf-dimensional queries in 0{n l ~ l l d + k) time. 



A number of geometric problems can be solved more efficiently when points lie on a grid, i.e., 
when coordinates of points are integers^] bounded by a parameter U. In the case of range reporting, 
significant speed-up can be achieved if the set of points S does not change. There are static data 
structures that support orthogonal range reporting queries in 0(loglogC7 + k) time [17|4] . On the 
other hand, if points can be inserted into or deleted from S, then any data structure that supports 
updates in log°^ n time needs i?(logn/ log logn + k) time to answer a two-dimensional range 
reporting query [5]. This bound is also valid in the case when all points belong to a U x U grid. 

Our Result. In this paper we consider the situation when coordinates of moving points belong 
to a Ux U grid. Our data structure supports orthogonal range reporting queries in O ( i/log U / log log U + 
k) time. This result is valid in the standard kinetic model with additional conditions that all points 
move with fixed velocities along linear trajectories and all changes in the trajectories are known in 
advance. Queries can be answered at any time t, where t is a positive integer bounded by U 0<k1 \ 
Updates are performed only when x- or y-coordinates of any two points in S swap their relative 
positions, and each update takes poly-logarithmic time. The total number of events after which the 
data structure must be updated is 0(n 2 ). For instance, for U = n ^ our data structure answers 
queries in O ( -\/log n / log log n + k) time. Our result also demonstrates that the lower bound for 
dynamic range reporting queries can be surpassed in the case when the set S consists of linearly 
moving points. Our data structure uses 0(n log 2 n) space and supports updates in 0(log 3 n) time, 
but space usage and update cost can be reduced if only special cases of reporting queries must be 
supported. We describe a 0(n) space data structure that supports updates in O(logn) time and 
dominance queries in 0(y / log U / log log U + k) time. We also describe a 0(n log n) space data struc- 
ture that supports updates in 0(log 2 n) time and three-sidec@ queries in 0{ yjlog U / log log U + k) 
time. 

2 Overview 

In section [3] we show that we can find the predecessor point of any v £ U in the set S (with respect 
to x- or y-coordinates) in 0{ y^log U / log log U) time by answering a point location query among 
a set of segments. In fact, identifying the predecessor of a point q is the bottleneck of our query 
answering procedure. 

In section H] we describe the data structure that reports all points p G S that dominate the 
query point q, i.e. all points p such that p.x > q.x and p.y > q.y; henceforth p.x and p.y denote the 
x- and y-coordinates of a point p. The query time of our data structure is 0(y / log U / log log U + k). 
The data structure is based on the modification of the d-approximate boundary [TO] for the kinetic 
framework. The d-approximate boundary [19] enables us to obtain an estimation for the number of 
points in S that dominate an arbitrary point q. If a point q dominates a point on a d-approximate 
boundary A4, then q is dominated by at most 2d points of S; if q is dominated by a point on A4, 
then q is dominated by at least d points of S. In section U] we show that a variant of a d- approximate 
boundary can be maintained under kinetic events. If a query point q is dominated by k < logn 
points of S, we can reduce the dominance query on S to a dominance query on a set that contains 
O(logra) points using a (/-approximate boundary for d = logn; see section|H Otherwise, if k > logn, 
we can answer a query in 0(logn + k) = O(k) time using a standard kinetic data structure [2]. 

For simplicity, we assume that all points have positive coordinates. 

The query range of a dominance query is a product of two half-open intervals. The query range of a three-sided 
query is a product of a closed interval and a half-open interval. 
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A data structure that supports dominance queries can be transformed into a data structure that 
supports arbitrary orthogonal range reporting queries |ll|18j by dividing the set S into subsets Si 
and constructing dominance data structures for each Si as described in section However, we may 
have to delete a point p from one subset Si and insert it into a subset Sj after a kinetic event. 
Unfortunately, the construction of |19| is static. It is not clear how (and whether) to modify the 
(i-approximate boundary, so that insertions and deletions are supported. However, in our case the 
deleted (inserted) point always has the maximal or minimal x- or y-coordinate among all points in 
Si (Sj). We will describe in section [5] how our dominance data structure can be modified to support 
these special update operations without increasing the query time. Our technique is similar to the 
logarithmic method and can be of independent interest. Thus we obtain the data structure for 
general orthogonal range reporting queries. 

3 One-Dimensional Searching 

Let S x (t) and S y (t) denote the sets of x- and y-coordinates of all points at time t. In this section 
we will describe how we can identify the predecessor of any q x in S x (t) (resp. of q y in S y (t)) at 
current time t in 0( y/log U / log log U) time using a linear space data structure. 

Let Xi(t) = a^t + bi be the equation that describes the x-coordinate of the point pi £ S at time 
t. The trajectory of the point in (t, x) plane (t-axis is horizontal) is a sequence of segments. Since 
we assume that all changes of point trajectories are known in advance, endpoints of all segments 
are known in advance. Two points swap ranks of their x-coordinates at time t if and only if their 
segments intersect at time t. We can find intersection points of all segments using the standard 
sweepline algorithm [10] in 0((n + /) logn) time, where / is the number of segment intersections. 
We start the sweepline at t = and move it to the right until n intersection points are identified 
or the last intersection point is found. These intersection points and the corresponding segments 
induce a subdivision of the (x, t) plane of size 0(n). We can construct the data structure for planar 
point location [12] that supports queries in 0(-y/log U / log log U) time. Let t max be the largest t- 
coordinate of the already processed intersection point. For t < t max , we can find the predecessor 
of any x by locating the segment lying immediately below the point (t,x). When t = t max , we 
continue the sweepline algorithm and find the next n segment intersection points. The algorithm 
described in [10J finds / next segment intersection points 0(/logn) time. Since the point location 
data structure [12] for a subdivision of size / can be constructed in O(f) time, an amortized 
cost of processing a kinetic event is O(logn). We can de-amortize the update cost using standard 
techniques. 

4 Dominance Queries 

In this section we describe the data structure that reports all points from S that dominate the 
query point q, i.e. all points in the region [q.x,+oo) x [q.y, +oo). Our data structure is based on 
maintaining the (i-approximate boundary for a set S. The notion of a d-approximate boundary is 
introduced in [19]; in this paper we change the definition and describe a kinetic version of this 
construction. 

For a horizontal segment s, we denote by start (s) and end(s) x-coordinates of the left and 
the right endpoint of s; we denote by y(s) the y-coordinate of all points on s. We will say that 
a segment s covers a point p if the x-coordinate of p belongs to [start(s), end(s)]. In this paper 



3 



Fig. 1. An example of a (i-approximate boundary for d = 6. 



we define a d-approximate boundary as a polyline A4 that consists of alternating horizontal and 
vertical segments, divides the plane into two parts, and satisfies the following properties: 

Invariant 1 Let s and r be two consecutive horizontal segments. Then \{p G S \ start(s) < p.x < 
end(r)}| > d/2. 

Invariant 2 When a new segment s is inserted, the left endpoint of s is dominated by at most 
3d/2 points of S. The number of points in S that dominate the left endpoint of a segment s G M 
does not exceed 2d. 

Invariant 3 When a new segment s is inserted, the right endpoint of s is dominated by at least d 
points. The number of points in S that dominate the right endpoint of a segment s G A4 remains 
constant. 

Using Invariants [HO we can prove the following Lemma. 

Lemma 1. Every point on a d- approximate boundary A4 is dominated by at least d points and at 
most 2d points of S. There are 0(n/d) horizontal segments in M. 

Proof: If a point p G Ai is dominated by k points from S, then the left endpoint of some segment 
s is dominated by at least k points and the right endpoint of some segment r is dominated by at 
most k points. Hence, it follows from Invariants [2] and [3] that d < k < 2d. By Invariant [IJ for two 
consecutive segments r and s there are more than d/2 points p G S, such that p.x belongs to the 
interval [start (r), end(s)]. For each point p, p.x belongs to at most two such intervals; hence, the 
total number of segments is less than 8n/d. □ 

An example of a d-approximate boundary is shown on Fig [TJ We will show below how the concept 
of a d-approximate boundary can be used to support dominance queries in 0( \/log U j log log U + k) 
time. Later in this section we will show how Invariants [T] -[3] can be maintained. 
Kinetic Boundary. We will use a kinetic variant of the d-approximate boundary, i.e. segments of 
the boundary move together with points of S. For every horizontal segment s in a boundary A4, 
let l(s) denote the point with the largest y-coordinate such that l(s).y < y(s) and let u(s) denote 
the point with the smallest y-coordinate such that u(s).y > y(s). Let first (s) denote the point 



with the smallest x-coordinate such that f irst(s).x > start(s); let last(s) denote the point with 
the largest x-coordinate such that last(s).x < end(s). We assume that start (s) = f irst(s) — \ 
and y(s) = u(s) — ^. That is, the left end and the y-coordinate of a segment change when u(s) 
and f irst(s) move. The right end of the previous segment and the x-coordinate of the connecting 
vertical segment change accordingly. 

Answering Queries. Our data structure is based on a d-approximate boundary At of S for 
d = log n. For each segment s G M we maintain the set Dom{s) of all points that dominate the left 
endpoint of s. Obviously, the set Dom(s) changes only when events concerning f irst(s) or u(s) 
take place. 

All points of Dom(s) are stored in a data structure D s that supports dominance queries in 
0(logd+k) = 0(loglogn + k) time. We can use the data structure of [2J, so that the space usage is 
0(d) and updates after events are supported in O(loglogn) time. It is possible to modify the data 
structure of [2J, so that points can be inserted into Dom(s) or deleted from Dom(s) in O(loglogn) 
time. All points of S are also stored in a kinetic data structure G that uses 0(n) space and supports 
dominance queries in 0(logn + k) time and updates after kinetic events in O(logn) time. Again, we 
can use the result of [2] to implement G. Finally, we must be able to identify for each point p E S 
the segment s 6 M. that covers p. Using the dynamic union-split-find data structure of [13] or the 
van Emde Boas data structure |14| . we can find the segment that covers any p 6 S in O (log log n) 
time. When a new segment is inserted into or deleted from A4, the data structure is updated in 
O (log log n) time. 

Given a query point q, we identify the point p with the largest x-coordinate such that p.x < q.x; 
this can be done in 0( y^og U / log log U) time using the construction described in section [3J The 
point q dominates a point on Jv[ if and only if q dominates the left endpoint of the segment s that 
covers p or the left endpoint of the segment h that follows s. Suppose that q dominates a point 
on A4. Then q is dominated by at most 2d points of S. Let v be the left endpoint of a segment 
g, such that v is dominated by q. Each point p € S that dominates q also dominates v; hence, all 
points that dominate q belong to Dom(g). We can use the data structure D g and report all points 
that dominate q in 0(log log ra + k) time. Suppose that q does not dominate any point on A4. Then 
there is at least one point on M that dominates q. Hence, there are k > d = logn points of S that 
dominate q. Using data structure G, we can report all points that dominate q in 0(log n + k) = 0(k) 
time. 

Maintaining the d- Approximate Boundary. It remains to show how to maintain Invariants [Ul3] 
after operations x-move and y-move. Henceforth, we will use the following notation. Suppose that 
p.x < q.x before some kinetic event and p.x > q.x after this event. Then, we say that p is x-moved 
behind q (q is x-moved before p). Suppose that p.y < q.y before some kinetic event and p.y > q.y 
after this event. Then, we say that p is y-moved above q (q is y-moved below p). Each kinetic event 
can be represented as a combination of at most one x-move and at most one y-move. 

First, we consider the Invariant[2j Suppose that the left endpoint of an interval s is dominated by 
2d points of S. Let h be the segment that precedes s, i.e., end(/i) = start (s). Let v be the vertical 
segment that connects h and s. We look for a point p with p.x = start(s) and y(s) < p.y < y(h) 
such that p is dominated by 3d/2 points. If such a point p exists, we set pi = p and search for a 
point p r with start (s) < p r .x < end(s) and p r .y = p\.y such that p r is dominated by d points of 
S. If there is no such point, i.e., if the point (end(s),pj.y) is dominated by at least d points, we 
replace s with a segment s' such that the left endpoint of s' is p\ and the right endpoint of s' is 
the point (end(s),p/.y); see Fig[2ji. In other words we change the y-coordinate of s to p.y. The new 
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Fig. 2. Updating the d-approximate boundary when the Invariant [2] is violated. For simplicity point 
of S are not shown, (a) There exists a point pi but there is no p r . (b) There exist both pi and p r . 
(c) There are neither pi nor p r . 



set Dom(s) contains 3d/2 points and can be constructed in 0(d) time. If there is a point p r with 
start(s) < p r .x < end(s) and p r .y = p\.y that is dominated by d points of S, then we replace 
s with two new segments s' and s" . The left and right endpoints of s' are p[ and p' r respectively. 
The left endpoint of s" is the point p" such that p".x = p r .x and p" .y = y(s). See Fig [2)3. The 
set Dom(s') contains 3d/2 points. There are at most d/2 points q such that q.x > start(s) and 
y(s) < q.y < p.y; hence, there are at most d/2 points q such that q.x > p r .x and y(s) < q.y < p.y. 
Therefore, since p r is dominated by d points of S, p" is dominated by at most 3d/2 points of S and 
Dom{s") contains at most 3d/2 points. Since p" is dominated by the right endpoint of the segment 
s, Dom{s") contains at least d points. We can construct Dom(s') and Dom(s") and data structures 
D s i and D s n in 0(d) time. 

If the right endpoint of h is dominated by at least 3d/2 points, we shift the vertical segment v 
in +x direction, so that the right endpoint of h is dominated by d points from S or the segment 
s is removed. That is, we identify the point r with r.y = y(h) such that either r.x = end(s) and 
r is dominated by at least d points of S or r.x < end(s) and r is dominated by d points of S. If 
r.x < end(s), we set end(h) = start(s) = r.x and update Dom(s), D s accordingly (0(d) points 
are removed from Dom(s) and D s ). See Fig[2k- The new left endpoint of s is dominated by at most 
3d/2 points. If r.x = end(s), we remove the segment s with D s and Dom(s). 

The update procedure removes at most one segment and inserts at most two new segments that 
satisfy the Invariant [2j Hence, the update procedure takes 0(d) = O(logn) time. 

Now we describe how the Invariant [T] can be maintained. Let h and s be two consecutive segments 
(end(/i) = start(s)) and suppose that there are d/2 points that belong to [start (h), end(s)]. We 
replace h and s with one new segment g as follows. If the left endpoint of h is dominated by at least 
3d/2 points, then we set end(/i) = end(s) and remove the segment s, Dom(s), and D s . The point 
q with q.y = y(h) and q.x = end(s) is the new right endpoint of h. Since there are at most d/2 
points p G S such that start (s) < p.x < q.x, q is dominated by at least d points of S. Hence, the 
new segment h satisfies Invariant [3j See Fig. [3k. If the point I with l.x = start(/i) and l.y = y(s) 
is dominated by at most 3d/2 points of S, we set start(s) = start(/i) and remove h, Dom(h) and 
Dfi. The set Dom(s) and the data structure D s are updated. See Fig. [3b. Since there are less than 
d/2 points p € S, such that p.y > y(s) and start(/i) < p.x < end(h), 0(d) new points are inserted 
into Dom(s) and D s . If the left endpoint of h is dominated by less than 3d/2 points and the point 
/ is dominated by more than 3d/2 points, then we replace h and s with a new segment g. The left 
endpoint of g is the point m such that m.x = start(/i), y(s) < m.y < y(h), and m is dominated by 



6 



h q h h 



s : s s 
1 j l ' | 1 

(a) (b) (c) 

Fig. 3. Updating the d-approximate boundary when the Invariant Q] is violated. For simplicity 
points of set S are not shown, (a) The left endpoint of h is dominated by > 3d/2 points, (b) The 
point I is dominated by at most 3d/2 points, (c) The left endpoint of h is dominated by less than 
3d/2 points but I is dominated by more than 3d/2 points. 

3d/2 points. The right endpoint of g is the point r with r.x = end(s) and r.y = m.y. See Fig. [3b. 
The point r is dominated by at least d points of S. Hence, g satisfies Invariants [2] and El 

Now we turn to Invariant [3j The number of points that dominate the right endpoint of a 
horizontal segment s changes: 1) if a point p with p.x > end(s) is y-moved above u(s) or below 
u(s); 2) if a point p with p.y > y(s) is x-moved before last(s) or behind last(s). 

First, we consider y-moves. Essentially, when a point p is y-moved, we shift the segment s in 
+y or — y direction so that the number of points that dominate the right endpoint of s remains 
unchanged. Suppose that a point p = u(s) with p.x > end(s) is y-moved below l(s); let y n denote 
the y-coordinate of p. 

(a) If start (s) < l(s).x < end(s), then we change the y-coordinate of s so that y(s) = y n — \. The 
old point l(s) is added to Dom(s) and D s . 

(b) If l(s).x < start (s), then we also change the y-coordinate of s so that 

y(s) =y n -\- 

(c) If l(s).x > end(s), then we change the y-coordinate of s so that y(s) = l(s).y — h. We delete p 
from Dom(s) and D s and insert l(s) into Dom(s) and D s . 

Suppose that a point p with p.x > end(s) is y-moved above u(s) 

(d) If start (s) < u(s).x < end(s), then we change the y-coordinate of s so that y(s) = q.y — 5, 
where q is the point with the smallest y-coordinate such that q.y > y n . The old point u(s) is 
removed from Dom(s) and D s . If q = u(/i) for the segment h that precedes s, then the segment s 
and Dom(s) are deleted. 

(e) If u(s).s < start(s), then we also change the y-coordinate of s so that y(s) = q.y — 5, where q 
is the point with the smallest y-coordinate such that q.y > y n . If q = u(/i) for the segment h that 
precedes s, then the segment s and Dom(s) are deleted. 

(f) If u(s).x > end(s), then we change the y-coordinate of s so that y(s) = y n — \- We delete u(s) 
from Dom(s) and D s and insert p into Dom(s) and D s . 

Observe that the number of points in D s remains unchanged or increases by 1. See Fig[H 

We can handle the x-moves in a similar way. Let h be the horizontal segment that follows s in 
Ai, i.e., end(s) = start(/i). When a point p is x-moved we change end(s) (and start(s)) so that 
the number of points that dominate the right endpoint of s remains unchanged. Suppose that a 
point p = last(s) with p.y > y(s) is x-moved behind f irst(/i); let x n denote the new x-coordinates 
of p and let q be the point with the smallest x-coordinate such that q.x > x n . 
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(c) 



(d) 



(e) 



(f) 



Fig. 4. The segment s is shifted so that Invariant [I] is maintained after a y-move. Figures (a), 
(b),(c), (d), (e), and (f) correspond to cases (a), (b), (c), (d), (e), and (f) respectively. 



(a) If y(h) < first (h).y < y(s), then we set start (/i) = end(s) = q.x — 5 and remove first (fo) 
from Dom(h) and D^. If q = f irst(Zi') for the horizontal segment h! that follows h, then we delete 
the segment h and Dom(h). 

(b) If f irst(h).y > y(s), then we set start(/t) = end(s) = x n — \\ we also remove f irst(/i) from 
Dom{h) and D/, and insert p into Dom(h) and D/j. 

(c) If first(/i).y < y(/i), then we set start(/i) = end(s) = g.x — ^. If q = first(Zi') for the 
horizontal segment h! that follows h, then we delete the segment h and Dom(h). 

Suppose that a point p with p.y > y(s) is x-moved before last(s). 

(d) If y(h) < last(s).y < y(s), then we set start(/i) = end(s) = x n — \\ we add last(s) to Dom(h) 
and Dfr. 

(e) If last(s).y > y(s), then we set start(/i) = end(s) = last(s).x — ^; we also remove p from 
Dom{h) and -D^ and insert last(s) into Dom(h) and 

(f) If last(s).y < y(h), then we set start (/i) = end(s) = x n — |. We add last(s) to Dom(h) and 
O h 

Again the number of points in D s remains unchanged or increases by 1. 

We will show in the full version that we update Mi because Invariants [2] or [1] are violated at 
most once for Q{d) events. Update procedures for maintaining Invariants Q] and [2] involve inserting 
and deleting a constant number of segments into A4 and the data structure D s for every segment s 
contains 0(d) points. Hence, we must perform 0(1) amortized updates of data structures D s after 
each kinetic event. Since every update of D s takes O(loglogn) time, the amortized cost of updating 
after an event is O (log log n). Since each data structure D s can be constructed in O(logn) time, 
the worst-case update time is O(logn). 

Construction of a d- Approximate Boundary. Now we show that d-approximate boundary can 
be constructed in 0(n) time if points are sorted by x- and ^-coordinates. Since a d-approximate 
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boundary consists of alternating horizontal and vertical segments, it suffices to determine the end- 
points of horizontal segments. We can guarantee that the left endpoint of each segment is dominated 
by 3d/2 points of S and the right endpoint of each segment is dominated by d points of S using 
the following algorithm. Lists L x and L y contain points of S sorted in descending order of their 
x-coordinates and y-coordinates respectively. With every element of L x we store a pointer to its 
position in L y and vice versa. The pointer ptr x (ptr y ) points to the first not yet processed element 
in L x {Ly). We construct a sequence of horizontal and vertical segments so that the left endpoint 
of each segment is dominated by 3d/ 2 points and the right endpoint of each segment is dominated 
by d points. 

We assign ptr y to the (3d/ 2 + l)-th element of L y and ptr x to the first element of L x ; the point 
Pi with p\.y = ptr y .y — 1/2 and p.l.x = is the left endpoint of the first segment. Clearly, p\ is 
dominated by 3d/2 points. (1) The right endpoint of the segment with left endpoint pi can be found 
as follows. We traverse elements of L x that follow ptr x until d/2 + 1 points q such that q.y > ptr^.y 
are visited and update ptr x accordingly. The point p r with p r .x = ptr^.x — 1/2 and p r .y = p\.y is the 
right endpoint of the currently constructed horizontal segment. (2) We identify the left endpoint 
of the next segment by traversing elements of L y that follow ptr^ until d/2 points q such that 
q.x > ptic x .x are visited or we reach the end of the list L y . The pointer ptr y is updated accordingly. 
If we reached the end of L y , then all points are processed and the algorithm is completed. Otherwise 
we set pi.x = ptr x .x — 1/2 and p\.y = pti y .y — 1/2, go to step (1) and determine the right endpoint 
of the next segment. 

Theorem 1. There exists a linear space data structure that supports dominance queries for a set 
of linearly moving points on U x U grid in 0(i/log U / log log U + k) time and updates after kinetic 
events in O(logn) worst-case time. The amortized cost of updates is O(loglogn). If trajectories of 
the points do not change, then the total number of kinetic events is bounded by 0(n 2 ). 

5 Orthogonal Range Reporting Queries 

Three-sided range reporting queries and orthogonal range reporting queries can be reduced to 
dominance queries with help of standard techniques. However to apply these techniques in our 
scenario, we must modify the data structure of section U so that insertions and deletions are 
supported in some special cases. At the end of this section we demonstrate how our data structure 
with additional operations can be used to support arbitrary orthogonal range reporting queries. 

Additional Update Operations. We first describe the data structure that supports insertions 
and deletions in two special cases: Let x m \ n and x max be the smallest and the largest x-coordinates 
of points in S. The operation insert^ inserts a point p with x max < p.x. The operation delete^ 
deletes a point p with p.x = x max . Operations insert" and delete" insert and delete a point whose 
^-coordinate is smaller than x-coordinates of all other points in S. It is easy to augment our data 
structure so that insert" and delete" are supported: the inserted (deleted) point is either below a 
(^-approximate boundary A4 or dominates only the leftmost horizontal segment of M. Hence, each 
insert" and delete" affects the data structure D s for at most one segment s. Essentially we can 
handle insert" and delete" in the same way as x-moves for the leftmost segment s. Maintaining 
the (i-approximate boundary M. after insert^" and delete^" is more involved: since the y-coordinate 
of a newly inserted (deleted) point can be larger than the y-coordinates of all (other) points in p, 
we may have to update data structures D s for all segments s € A4 after a single update operation. 
Below we describe how insert^ and delete^" can be supported. 
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Our approach is similar to the logarithmic method |9|16| that is used to transform static data 
structures into data structures that support insertions. We construct the data structure of section 0] 
augmented with insert" and delete" for sets H2,H%, . . . ,H m . A point p x-overlaps with point q if 
p.x > q.x. A point p x-overlaps with a set S if it x-overlaps with at least one point q G S. Each set 
H{ satisfies the following conditions: 

1. For i = 2, . . . ,m — 1, Hi contains between 2 2 * -1 and 2 2j+2 points; H m contains between 2 2 ' -2 
and 2 2i+2 + 2 2i points 

2. Each point of Hi x-overlaps at most 2 2l_4 points in Hi_\ 

3. At most 2 2l_3 points from Hi x-overlap with 

As follows from conditions 2 and 3, no element of Hi + \ x-overlaps with Hi-\\ The rightmost point 
in Hi + \ x-overlaps at most 2 2 * -2 leftmost points in Hi by condition 2. Only 2 2 * -3 rightmost points 
in Hi can x-overlap with a point in Hi-\. Hence, any q € Hi + \ that x-overlaps a point in 
would x-overlap \Hi\ — 2 2i ~ 3 > 2 2i ~ 2 points in Hi, which contradicts condition 2. 

For j = 1, ... ,m, we maintain max^j) = m&x{p.y\p € Hj} and mm x (j) = mm{p.x\p € Hj}. 
All minj;(j) are stored in a kinetic binary search tree. For each j, 2 < j < m, a point pj such that 
Pj.y = maxy(j) and pj.x = j is stored in a kinetic data structure Y; since Y contains O(logre) 
points, Y supports dominance reporting queries in 0(loglogn + k) time. The data structure Dj, 
j = 2, ... ,m, contains all points from Hj and supports dominance queries, x-moves and y-moves as 
described in section HI To speed-up update operations, we store only one data structure G for all 
points in \J^ =2 Hj. G is implemented as described in [2j, supports dominance queries in 0(logn + k) 
time and can be modified to support arbitrary updates as well as kinetic events in O(logra) time. 
We also store one data structure V that contains x-coordinates of all points If^^Hj and enables us 
to search in the set of x-coordinates at any time d. That is, all Dj share one data structure G and 
one data structure V. For each j we also store all points of Hj in a list Lj that contains all points 
from Hj in the descending order of their y-coordinates. 

Given a query Q = [a, +oo) x [b, +oo), we can find in O(loglogn) time the smallest index j, 
such that at least one point in Hj has x-coordinate smaller than a. Then, as follows from conditions 
2 and 3, x-coordinates of all points in H2 U . . . U i2j-i are greater than a, and both Hj and Hj + \ 
may contain points whose x-coordinates are greater than or equal to a. Sets Hf, f > j + 1, contain 
only points whose x-coordinates are smaller than a. Using Y, we can identify all Hf such that f < j 
and Hf contains at least one point p with p.y > b. For every such / all points p such that p € Hf 
and p.y > b can be reported by traversing the list Lf. Hence, reporting all points p G Hf such that 
/ < j, p.x > a and p.y > b takes 0(loglogn + k) time. We can report all points in Hj and Hj + \ 
that belong to Q = [a, +00) x [b, +00) in O^log U/ log log U + k) time using data structures Dj 
and Dj + \ respectively. 

It remains to show how conditions 1-3 above can be maintained. Clearly, conditions 1-3 are 
influenced by x-moves and operations insert" , delete" , insert^" and delete^" ; y-moves cannot violate 
them. We say that a set S is x-split into sets S\ and S2 if S\ U S2 = S and the x-coordinates of 
all points in S\ are larger than the x-coordinates of all points in S2. After an operation insert^ or 
delete^", we re-build the data structure for H2 in O(l) time. When the number of elements in a set 
Hj, j < m, becomes smaller than 2 2 - 7-1 or greater than 2 2j+1 , we x-split the set Hj + \ U Hj into two 
new sets H'- and H'- +l , so that H'- contains 2 2j points and H'- +l contains + — 2 2j points. 

If the number of elements in H m exceeds 2 2m+2 + 2 2m we x-split Hj into sets H' m and H' m+1 that 
contain 2 2m and 2 2m+2 points respectively. Suppose that the number of points in H m is smaller 
than 2 2m ~ 2 . If |J5T ro _i| + \H rn \ < 3 • 2 2m " 2 , we set H' m _ x = H m _ x U H m and decrement m by 1. 
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If \H m _i \ + \H m \ > 3 • 2 2m " 2 , we x-split H m - X U H m into H' m _ x and H' m that contain 2 2m " 2 and 
|ff m _i| + |i/ m | — 2 2m ~ 2 points respectively. We also take care that conditions 2 and 3 are maintained. 
If, as a result of x-moves, the number of points in some Hj that x-overlap Hj—\ exceeds 2 2j ~ 3 , or 
there is at least one point in Hj that x-overlaps more than 2 2j ~ 4 points in iij-i, then we x-split 
U Hj into sets and H'- that contain and \Hj\ elements respectively. Each set Hj 

is rebuilt at most once after a sequence of &{Hj) special insert or delete operations. Each Hj is 
also re-built at most once after a sequence of 0(\Hj\) x-moves, i.e., after 0{Hj) kinetic events. We 
can maintain the list of points in each Hj sorted by their y-coordinates; hence, the data structure 
Dj for a newly re-built set Hj can be constructed in 0(\Hj\) time. Therefore, the amortized cost of 
updates and kinetic events is O(logn). We can de-amortize update costs using the same techniques 
as in the logarithmic method [16J. 

Lemma 2. There exists a 0(n) space data structure that supports dominance queries on U x U 
grid in 0{ y^og Uf log log U + k) time. Updates after kinetic events and operations insert^ - ; delete^", 
insert" and delete" are supported in O(logn) time. 

Three-Sided Reporting Queries. Now we are ready to describe data structures that support 
three-sided reporting queries, i.e., the query range is a product of a closed interval and a half- 
open interval. We apply the standard method used in e.g. [19], [18] to augment data structure for 
dominance queries. 

Let T be an arbitrary balanced tree with constant node degree on the set of x-coordinates of 
all points in S. We associate an interval (ai,bj\ with each leaf I of T, where a/ is the predecessor 
of the smallest value mi stored in the node I, and bi is the largest value stored in the node I. 
We associate an interval (a;,+oo) with the rightmost leaf I. With each internal node v of T we 
associate an interval int(u) = Uint(t> j) for all children Vi of v. Let S v be the set of points p £ S such 
that p.x € int(u). In every internal node v we store two data structures C v and 1Z V that support 
dominance reporting queries open to the left and open to the right (i.e., queries (— oo,a] x (—00,6] 
and [a, +00) x (—00, b\) for the set S v . In each node v we also store the list of points in S v sorted 
by their y-coordinate. 

Given a three-sided query with Q = [a, b] x (—00, c] we can find in time 0(y / log U / log log U) the 
node v such that [a, b] C int(u), but [a, b] (jL int(uj) for all children v% of v. Suppose ±nt(vj) C [a, b] 
for j = r, r, +1, . . . , q. Then x-coordinates of all points in children v r , . . . , v q of v belong to [a, b]. 
We can report all points p in v r , . . . , v q whose y-coordinate do not exceed c using sorted lists 
of points in S Vr , . . . , S Vq . We also answer two dominance queries Q\ = [a, +00) x (— 00, c] and 
Q2 = (—00,6] x (— 00, c] with help of data structures 7£„ r _i and C Vq+1 respectively. 

After a kinetic event affecting two points p and q from the same set S v , the data structures C v 
and Tl v are updated. After a kinetic event that affects points p and q that belong to two neighbor 
sets S Vi and S Vi+1 respectively, we swap p and q: p is removed from S Vi and inserted into S Vi+1 , 
and q is removed from S Vi+1 and inserted into S Vi . In this case a constant number of operations 
insert" and delete" (resp. insert^" and delete^") is performed. Each point belongs to O(logn) sets 
S v . Hence, the space usage is 0(n log n) and an update after a kinetic event takes 0(log 2 n) time. 

Lemma 3. There exists a 0(n log n) space data structure that supports three-sided reporting queries 
for a set of linearly moving points on a U x U grid in 0(y / log XJj log log U + k) time and updates 
after kinetic events in 0(log 2 n) time. If trajectories of the points do not change, then the total 
number of kinetic events is 0(n 2 ). 
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Orthogonal Range Reporting Queries. In a similar way to Lemma [2] we can extend the data 
structure to support update operations in two other special cases. Let y m i n and y max be the smallest 
and the largest y-coordinates of points in S. The operation insert,}" inserts a point p with y max < P-V- 
The operation delete^" deletes a point p with p.y = y max . Operations insert" and delete" insert and 
delete a point whose y-coordinate is smaller than y-coordinates of all other points in S. Again, it 
is easy to modify the data structure of Lemma [3] so that it supports insert" and delete" because 
these operations affect at most one segment of Ai. We can support insert^ and delete^" using the 
same construction as in Lemma El 

A point p y-overlaps with a point q if p.y > q.y. Analogously to Lemma O points are stored in 
sets V2, • • • , V m and we maintain the invariants: 

1. For i = 2, ... ,m — 1, Vi contains between 2 2 *" 1 and 2 2 * +2 points; V m contains between 2 2i ~ 2 
and 2 2i+2 + 2 2i points 

2. Each point of Vi y-overlaps at most 2 2 *" 4 points from V%-\ 

3. At most 2 2 *" 3 points from Vi y-overlap with Vi-i 

Elements of Vi are stored in an augmented data structure Ei of Lemma [2] so that kinetic 
events and operations insert", delete", insert^, delete^, insert" and delete" are supported. For 
j = 1, . . . , m, we maintain max x (j) = max{p.x\p £ Vj} and min^j) = mm{p.y\p € Vj}. All mm y (j) 
are stored in a kinetic binary tree. For each j, we store a point pj with pj.x = max z (j) and pj.y = j 
in a kinetic data structure X; X supports dominance reporting queries in 0(loglogn + k) time 
and updates in O(logra) time. For each j we also store all points of Vj in a list L'j that contains all 
points from Vj in the descending order of their x-coordinates. 

Given a query Q = [a, +00) x [b, +00), we can find in O(loglogn) time the smallest index j, 
such that at least one point in Vj has y-coordinate smaller than b. According to conditions 2 and 
3 above, Vj and Vj+\ may contain points whose y-coordinate are greater than or equal to b. The 
y-coordinates of all points in V2 U ... U Vj—\ are greater than b. The y-coordinates of all points in 
sets Vi, i > j + 1, are smaller than b. 

Using X, we can identify all Vf such that f < j and Vf contains at least one point p with p.x > a. 
For every such /, all points p such that p G Hj and p.x > a can be reported by traversing the list L'^. 
Hence, reporting all points p £ V/ such that f < j, p.x > a and p.y > b takes 0(log log n + k) time. 
We can report all points in Hj and Hj + \ that belong to Q = [a, +00) x [b, +00) in 0(log log U + k) 
time using data structures Ej and Ej + \ respectively. 

When a point p is inserted with an operation insert^ or insert", we identify the data structure 
Ei, such that p.y > min y (i) but p.y < min^j) for all j < i, and insert p into Ei as described 
in Lemma [2j When a point p is deleted with operations delete^ or delete", we delete p from the 
data structure Ei such that p £ Hi. When a point is inserted or deleted with operations insert" 
or delete", we insert or delete this point into the set V m . If a point is inserted or deleted with 
operations insert^" or delete^", we rebuild the data structure E%. Invariants 1-3 for sets Vi can be 
maintained under insert^, delete^, and kinetic events with the same method that was used in 
Lemma [2] to maintain sets Hi. 

Hence, we obtain 

Lemma 4. There exists a data structure that supports dominance queries on U x U grid in 
0(-^/log U / log log U + k) time. Kinetic events and operations insert^, delete^", insert", delete", 
insert^", delete^ , insert", and delete" are supported in O(logn) time. 
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We can obtain a 0(n log n) space data structure for three-sided queries that supports operations 
insert^" , delete^", insert", and delete" as well as kinetic events in 0(log 2 n) time in the same way 
as in the proof of Lemma [3j Using the same technique once again, we obtain the result for general 
two-dimensional range reporting queries stated in Theorem [2 

Theorem 2. There exists a 0(n log 2 n) space data structure that supports orthogonal range re- 
porting queries on U x U grid in 0(y / log U / log log U + k) time and updates after kinetic events in 
0(log 3 n) time. If trajectories of the points do not change, then the total number of kinetic events 
is 0(n 2 ). 

6 Conclusion and Open Problems 

In this paper we describe a data structure for a set of moving points that answers orthogonal range 
reporting queries in 0{^/\og U / log log U +k) time. The query time is dominated by the time needed 
to answer a point location query; the rest of the query procedure takes 0(log log n + k) time. Thus 
a better algorithm for the point location problem would lead to a better query time of our data 
structure. Proving any ,!?(log log [/) lower bound for our problem, i.e., proving that kinetic range 
reporting is slower than static range reporting, would be very interesting. 

While kinetic data structures usually support arbitrary changes of trajectories and only re- 
quire that points follow constant-degree algebraic trajectories, we assume that points move linearly 
and all changes of trajectories are known in advance. An interesting open question is whether we 
can construct a kinetic data structure that achieves o(logn/loglogre) query time without these 
additional assumptions. 
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